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Abstract. A substantial part of the dark matter of the Universe could be in the form of compact objects (MACHOs), detectable 
through gravitational microlensing effects as they pass through the line of sight to background light sources. So far, most at- 
tempts to model the effects of high-redshift microlensing by a cosmologically distributed population of MACHOs have assumed 
the compact objects to be randomly and uniformly distributed along the line of sight. Here, we present a more realistic model, in 
which the MACHOs are assumed to follow the spatial clustering of cold dark matter. Because of sightline-to-sightline variations 
in surface mass density, this scenario leads to substantial scatter in MACHO optical depths, which we quantify as a function 
of source redshift. We find that while optical depth estimates based on a uniform line-of-sight distribution are reasonable for 
the highest-redshift light sources, such estimates can be incorrect by a factor of «; 2 for the nearby (z ss 0.25) Universe. Hence, 
attempts to derive the cosmological density of MACHOs from microlensing observations of only a few independent sight- 
lines can be subject to substantial uncertainties. We also apply this model to the prediction of microlensing-induced variability 
in quasars not subject to macrolensing, and demonstrate that relaxing the assumption of randomly and uniformly distributed 
MACHOs only has a modest impact on the predicted light curve amplitudes. This implies that the previously reported problems 
with microlensing as the dominant mechanism for the observed long-term optical variability of quasars cannot be solved by 
taking the large-scale clustering of dark matter into account. 
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' 1. Introduction Griest [199T1 l, several non-baryonic dark matter candidates^ can 



also manifest themselves in this way - e.g. axion aggregates 



X 



Q | The identification of dark matter, estimated to contribute (Membrado 1998), mirror matter objects (Mohapatra& Teplitz 

around 90% to the total matter content of the Universe, re- 119991) . primordial black holes (Green 120001 1. quark nuggets 

mains one of the most important quests of modern cosmol- (Chandra & Goval l2000"b . preon stars (Hansson & Sandin 2005) 

ogy. To complicate the problem, dark matter appears to ex- and scalar dark matter miniclusters (Zurek et al. 120071 . Hence. 

\ ist in at least two different forms: one baryonic and on e non- MACHOs could in principle make up a substantial fraction of 

baryonic. According to current inventories (Fukugita 120041 the dark matter of the Universe. 

Fukugita & PeeblesEBa. ~ 1 /3 of the baryons (Q M ,baryo„ic * Potential MACHOs have indeed been detected through mi- 

0.04) are still unaccounted for, whereas nearly all (save for a crolen sing effects towards stars in the LMC (Alcock et al. 

small contribution from neutrinos) of the non-baryonic matter [200i3h . M31 (Calchi Novati et al. 120051) and towards multiply- 

(£Vnon-baryonic = q m - ^M.baiyonic ~ 0.2; e.g. Spergel et al. imaged quasars (e.g. Pelt et al. 119981 Koopmans & de Bruyn 

HHD remains elusive. 120001 Colley & Schild|2003). Claims of MACHO microlens- 

MACHOs (Massive Astrophysical Compact Halo Objects) ing of ga mma-ray bursts (Garnavich et al. 120001 Baltz & 

represent one class of dark matter candidates, which may be Hui |200 5J, and quasars whic h are n ot multiply-imaged (e.g. 

detected through gravitational microlensing effects as they pass Hawkins[l996l Zakharov et al. 120041) have also been made, but 

through the line of sight to distant light sources (e.g. stars, su- in these cases, the interpretation is more ambiguous, since in- 

pernovae, gamma-ray bursts, quasars). Although the MACHO trinsic variability cannot easily be ruled out. 
acronym was originally invented with baryonic objects like 

faint stars and stellar remnants in mind (Carr & Primack 1990; i T ■ " , , . t t r 

v 1 ' 1 In this category, one usually also includes candidates that are for- 
mally baryonic (like quark nuggets), but which evade the current con- 
straints on the cosmic baryon density and are expected to behave like 
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For the detections of MACHOs towards the LMC, alterna- 
tive explanations like self-lensing (i.e. microlensing by normal 
stars inside the target galaxy, rather than MACHOs in the dark 
halo of the Milky Way), variable stars or background super- 
novae have been suggested (e.g. Evans & Belokurov 2004]). It 
is also possible that blending effects introduced by the use of 
faint source stars may have affected the results in some unde- 
sired way (Tisserand et al. 2007 ). While these explanations can- 
not account for the microlensing events seen towards multiply- 
imaged quasars, the MACHO interpretation of some of the lat- 
ter observations has nonetheless been questioned, as the poorly 
understood structure of these light sources make it very difficult 
to robustly determine the nature of the microlenses in this lens- 
ing situation (e.g. Wyithe & Loeb 2002; Lewis & Gil-Merino 
[20061 . 

Although the local microlensing surveys impose strong up- 
per limits on the dark matter contribution from MACHOs in 
the 10~ 7 Mq-10 1 M q mass range (e.g. Lasserre et al. 2000; 
Tisserand et al. 120071 1. these constraints apply only to classes 
of MACHOs to which these surveys are sensitive. In fact, lo- 
cal microlensing observations are only able to detect the dens- 
est, most compact MACHOs at a given mass (e.g. Walker & 
Wardle [19991 Zackrisson [20051 Zurek et al. 120071 whereas 
more diffuse objects will give rise to detectable microlensing 
effects only at higher redshifts. 

Quasar emission-line and variability statistics also impose 
a number of important constraints on the cosmological density 
of MACHOs (Schneider [T993l Dalcanton et al. [19941 Wiegert 
2003 ), but these are sensitive to assumptions concerning the 
typical size of the optical continuum-emitting region in quasars 
and can in principle be sidestepped if this parameter is taken to 
be sufficiently large (Zackrisson & Bergvall 2003 ; Zackrisson 
et al. I2004I I. Disturbingly, the recent analysis by Pooley et al. 
(2006 2007 ) of combined X-ray and optical data on multiply- 
imaged quasars suggests that their optical continuum-emitting 
regions may be much larger than is usually assumed. 

Thus, the existence and relative importance of MACHOs 
is still a matter of debate. To remedy this unsatisfactory situa- 
tion, a large number of new microlensing surveys at both low 
and high redshift are currently being planned or are already 
underway. While the SuperMACHO project is continuing to 
monitor the LMC, the Angstrom, MEGA, POINT- AGAPE and 
WeCapp collaborations are targeting M31. Upcoming high- 
resolution facilities like the SIM PlanetQuest satellite also 
show great promise for the detection of astrometric microlens- 
ing effects, which will help to constrain the properties of the lo- 
cal MACHOs (e.g. Paczynski[[998] Han & Kim[l999] Rahvar 
& GhassemiEOQl]). 

To assess the possibility of more diffuse MACHOs and 
MACHOs on larger (galaxy-cluster) scales, microlensing mon- 
itoring of M87 in the Virgo cluster (Baltz et al. 2004) has 
started, and several other high-redshift MACHO searches 
based on cluster-quasar microlensing (Tadros et al. 2001) and 
cluster-cluster microlensing (Totani 2003) have also begun. 
With future supernova satellites like SNAP, JEDI or Destiny, 
it should also be possible to impose strong constraints on 
MACHOs from effects on the distribution of type la super- 
nova peak luminosities (e.g. Minty et al. 120021 Metcalf & 



Silk 120071 ) and by the temporal light curve fluctuations (e.g. 
Bagherpour et al. 2006 ; Dobler & Keeton 2006 ) that microlens- 
ing by foreground objects is expected to produce. Microlensing 
by MACHOs can also leave an imprint in the afterglows of 
gamma-ray bursts (e.g. Wyithe & Turner 2002b; Baltz & Hui 
120051 ). 

The interpretation of microlensing effects at cosmological 
distances is much more challenging than for local light sources. 
Due to the higher optical depth, there is a non-negligible 
probability that several microlenses may simultaneously con- 
tribute to the magnification. In the case of multiply-imaged 
light sources, it is customary to adopt the thin lens approxi- 
mation, as most of the microlensing can be assumed to take 
place in the dark halo responsible for the macrolensing ef- 
fect. For high-redshift light sources (i.e. quasars, supernovae, 
gamma-ray bursts) which are not multiply-imaged, the situa- 
tion becomes more complicated, as compact objects at vastly 
different distances along the line of sight may contribute to the 
microlensing. Most studies of the latter situation (e.g. Rauch 
[T99T1 Schneider[l993j Loeb & PernaHMl Minty et al. l2002l 
Zackrisson & BergvaU l2003l Zakharov et al. 120041 Baltz & Hui 
2005 ) adopt the Press & Gunn ( 119731 ) approximation, in which 
the microlenses are assumed to be randomly and uniformly dis- 
tributed with constant comoving density along the line of sight. 
Due to the strong clustering of matter, this approximation may 
not be entirely suitable and may give rise to misleading results 
(Wyithe & Turner 2002a). In this paper, we relax this approxi- 
mation in favour of a model where the compact objects instead 
follow the spatial distribution of cold dark matter, thereby clus- 
tering into cosmologically distributed halos and subhalos, as 
would be expected for non-baryonic MACHOs. We demon- 
strate that this gives rise to substantial sightline-to-sightline 
scatter in microlensing optical depths, and quantify this scat- 
ter as a function of source redshift. By applying the model to 
the microlensing contribution to the long-term optical variabil- 
ity of quasars which are not multiply-imaged, we also show 
that relaxing the Press & Gunn approximation only has a mod- 
est effect on the prediction of light curve amplitudes. Hence 
the conclusion reached by Zackrisson et al. (2003), who ar- 
gued that the long-term optical variability of quasars could not 
primarily be caused by microlensing - contrary to claims by 
Hawkins (e.g. [19961 [20001 [200T1 120021 ) - appears to be robust 
to the line-of-sight variations in matter content predicted in the 
ACDM cosmology. 

The matter distribution model and its numerical implemen- 
tation are described in Sect 2. In Sect. 3, we apply this model 
to the calculation of the microlensing optical depth for high- 
redshift light sources, and explore the impact of the various 
components of the model. In Sect. 4, we couple the matter 
distribution model to the microlensing code by Zackrisson & 
Bergvall (2003 ), and apply it to problem of predicting the mi- 
crolensing contribution to the long-term optical variability of 
quasars. Sect. 5 discusses the simplifications adopted in this 
work and a number of ideas for further refinement. Our find- 
ings are summarized in Sect. 6. 
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2. The distribution of dark matter along the line of 
sight to light sources at cosmological distances 

The matter distribution model used here is based on the ACDM 
scenario, where we assume Qm = 0.3, Qbar = 0.04 £2a = 0.7, 
cr 8 = 0.9 and H = I00h km s~' Mpc" 1 with h = 0.72. 
The overall results, however, are not sensitive to minor adjust- 
ments of these parameters, like the ones suggested by the 3- 
year WMAP data (Spergel et al. 120071 . 

2.1. MACHOs 



As the concept of MACHOs (Carr & Primack [19901 Griest 
I199U was introduced before the emergence of the ACDM 
model, and before non-baryonic matter had reached its current 
stature, it is not obvious how these objects should be defined, 
or how they fit into the hierarchical picture of modern cosmol- 
ogy. Should only point-mass objects, or objects with densities 
above a certain absolute threshold be considered MACHOs? 
Do MACHOs populate all halos, or just halos of a certain mass 
(e.g. galaxy-mass halos)? What about the low-density material 
that cannot easily be identified as belonging to any halo? 

Throughout this paper, we use the MACHO acronym to 
describe any kind of dark matter objects sufficiently compact 
to cause microlensing effects at cosmological distances, even 
though this definition converts into a distance-dependent den- 
sity threshold (e.g. Zackrisson I2005l l. We furthermore assume 
that MACHOs follow the spatial distribution of CDM, even 
when some fraction of the CDM is assumed to be smoothly 
distributed, i.e. either not associated with halos, or associated 
with halos below the mass resolution limit of our numerical 
model. 

2.2. Dark halos 

We assume that a redshift-dependent mass fraction /haio(z) 
of all matter is associated with dark matter halos of masses 
M in the interval M m i n < M < M max , whereas a fraction 
/smooth(z) = 1 - /haio(z) is smoothly distributed. We further- 
more assume that MACHOs constitute a scale- and redshift- 
independent mass fraction /macho = ^macho/^m of the mat- 
ter in both of these components, as the smooth component may 
indeed contain halos at masses below the adopted resolution 
limit. From now on, we reserve the term halo for objects that 
are not located inside other halos, and the term subhalos for 
those that are. For the mass function and spatial distribution 
of the halos, two options are implemented. Either a halo cata- 
logue from an N-body simulation is used, providing both the 
masses and positions of the halos, or the halos are assumed 
to be randomly distributed throughout the simulated volume, 
obeying the redshift-dependent mass function given by Sheth 
& Tormen ( [19991 : 



dn(M, z) 



p(z) diner 1 12a 

~ /i \ j 



d\ogM M dlogM V n 



5 

— exp 



-aS 2 



2cr 2 



(1) 



where n(M, z) is the number density of halos with masses below 
M at redshift z, p(z) is the mean matter density of the Universe, 
A = 0.3222, p = 0.3 and a = 0.707. 



Each dark matter halo is assumed to have a smooth, spher- 
ical component obeying the Navarro et al. ( 1 19961 119971 here- 
after NFW) density profile, given by: 



P(r) = 



(r/r s )(l+r/r s )2' 



(2) 



where rs is the characteristic scale radius of the halo and p\ is 
related to the density of the Universe at the time of collapse. 
The halos are truncated at the virial radius r v j r (z), defined as 
the radius at which mean enclosed density equals A v ; r z times 
the critical density p zz of the Universe at redshift z: 



7Vir(z) = 



3M 



1/3 



(3) 



\4n A vil -,, p c 

where A v ; r z can be approximated by (Bryan & Norman [1998l l: 
A vir ,„ « 18tt 2 + 82x- 39x 2 , (4) 
with x = O m (z) - 1 . 

The virial radius is linked to rs by the concentration param- 
eter c = r v i r /rs, with a median value of (Bullock et al. 2001 ): 



Mh 



1 + Z * iu - 1V1Q j 

and a scatter given by the log-normal distribution: 



pic) 



1 



V2^ 



exp 



C cr lnc 



(In c — In c) 2 

Inc 



(5) 



(6) 



where the standard deviation is cn nc = 0.3 (Wechsler et al. 
120021 . 

Given c, the pi parameter can be computed, using: 

_ Kir,z Pc,z 

Pi 3 ln(l +c)-c/(l +c)' 



2.3. Dark halo substructure 

Within each dark matter halo, we assume that a fraction 
/subhaiofc M) of the halo mass M is in the form of subhalos. To 
ensure mass conservation in the parent halo, the density p(r) of 
its smooth NFW component is decreased throughout the halo 
by a factor (1 - / su bhaio). For the subhalo mass distribution, we 
adopt the semianalytical model developed by van den Bosch et 
al. (120051) . in which the average mass fraction in subhalos for a 
parent halo of mass M is: 



log /subhalo = V0-41og(M/M*) + 5-2.74, 



(8) 



Here, M* represents the redshift-dependent characteristic mass 
scale defined by cr(M*, z) = S c (z), where cr 2 (M+, z) is the mass 
variance of the smoothed density field (which is related to the 
CDM power spectrum; see e.g. Eisenstein & Hu ll 9991 1 and <5 C 
is the critical threshold for spherical collapse (e.g. Navarro et 
al. [T9971 . 

In the van den Bosch model, the subhalo mass function for 
a parent halo of mass M is a Schechter function of the form: 



dN y (W\ I ¥ 

dln^F ~ pT(l -a) \~j3) 6XP \ /3 



(9) 
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where *P = M sub /M, p = 0.13, a = 0.996 - 0.028 log(M/M*) 
Here, y is given by: 

./subhalo 



y 



(10) 



P(l - a, IIP) - P(l - a, lO" 4 //?)' 

and P represents the incomplete Gamma function. The subhalo 
mass function <j9j is assumed to be valid throughout the mass 
range 10~ 4 < ¥ < / su bhaio- 

Due to the strong tidal stripping experienced by subhalos 
and the very high resolution required by N-body simulations 
to realistically treat this effect, the most appropriate form for 
the typical density profile of a subhalos is somewhat unclear 
at the moment. In the gravitational lensing literature, subhalos 
are usually treated as either NFW spheres (e.g. Metcalf 2005 ), 
singular isothermal spheres (e.g. Keeton 2003) or pseudo-Jaffe 
spheres (e.g. Dalai & Kochanek 2002). For simplicity, we as- 
sume here that the NFW profile adequately describes both ha- 
los and subhalos. As we will argue that subhalos represent 
one of the least important ingredients in the matter distribution 
model, this choice is not critical. 

Once the number and masses of the subhalos inhabiting a 
given parent halo have been determined using (0 and (0, these 
objects are distributed within the parent halo with a radial dis- 
tribution of the type suggested by Gao et al. (2004): 

(1 +acy 



n(< r) 



(11) 



N 1 + acr a ' 
where r is the distance to the parent halo centre in units of r v i r , 
n(< r) is the number of subhalos within r, N is the total number 
of subhalos within r v ; r , c is the concentration parameter of the 
parent halo, a = 0.244, a = 2 and p = 2.75. Only one gen- 
eration of subhalos is assumed, i.e. the possibility of subhalos 
within subhalos (e.g. Shaw et al. 120071) is neglected. 

2.4. Computational method 

To predict the matter distribution statistics for an ensemble of 
independent sightlines to the high-redshift Universe, a large 
number of sightline simulations are generated in a Monte Carlo 
fashion. To model the spatial distribution of dark matter be- 
tween the observer and source, this volume is divided into i 
lens planes, equally spaced in redshift. For each lens plane, a 
cylindric volume with radius R cy \ is considered, which is taken 
to be large enough to encompass the relevant features of the 
large-scale structure of the Universe. In practice, we take this 
cylinder to be equal to the virial radius of the most massive halo 
in the local Universe: R cy i = ^ v ii-,z=o(-^max)- Numerical tests in- 
dicate that choosing a simulation radius larger than this does 
not affect the statistical results. The total mass associated with 
each such cylinder is i taken to be: 

M cyU = Q M (1 + zd 3 p c 7rR 2 cyl AD<; (12) 

where p c is the critical density of the Universe and AD C is the 
light-travel time distance between the lens planes i - i and i. In 
a flat Friedmann-Robertson- Walker Universe (Q« + Q.\ =1), 
AD 1 : is given by: 

dz 



For each lens plane i, a fraction /haio,i of the mass M cy i ,- is as- 
sumed to be locked up in dark matter halos, whereas a fraction 
/smooth,/ = 1 - /halo,/ is taken to be smoothly distributed through- 
out the cylinder segment i. 

In scenarios involving halos with random spatial distribu- 
tion, each cylindrical volume i is populated with halos by ran- 
dom sampling of the dark halo mass function described in Sect. 
2.2 until the total halo mass reaches /haio./A/cyi,;- 

In scenarios involving spatially clustered halos, each cylin- 
drical volume i is inserted with random impact parameter into 
the N-body simulation cube closest in redshift, and the proper- 
ties (masses and coordinates) of the halos that happen to be lo- 
cated within that volume extracted. In cases where the N-body 
cube turns out to be smaller than the length D c of the cylinder, 
the part of the cylinder falling outside the cube is inserted into a 
new cube with a new random impact parameter. This procedure 
is repeated until every cylinder segment is completely located 
inside an N-body simulation cube. 

Once the masses M and positions of each halo inside each 
cylinder segment have been determined, a fraction / sub haio(^ Z) 
of each halo mass is redistributed in the form of subhalos 
through random sampling of the subhalo mass function and ra- 
dial distribution described in Sect. 2.3. 

Finally, we integrate the mass density along an infinitesi- 
mal beam through each cylinder centre, to derive the surface 
mass density of matter in that plane compared to the cosmic 
average E, at that redshift z,: 



f - 2 '- 



S/ 



% Q M (1 +Zi) 3 p ADf 



(14) 



AD' 



Ho i 



(l+z) V«m(1+z) + ^ 



(13) 



The stochastic variation of this quantity along each sight- 
line and between different sightlines is then used to predict 
sightline-to-sightline scatter in microlensing properties. 

2.5. Matter distribution scenarios 

To assess the importance of the various ingredients of the mass 
distribution model, five different levels (A-F) of refinement are 
tested. 

- A. All matter is assumed to be randomly distributed with 
constant comoving density along the sightline, i.e. accord- 
ing to the Press & Gunn ( 1973 ) approximation. 

- B. All matter is assumed to be in the form of smooth NFW 
halos (/halo = 1, /ub = 0) of mass 10 10 < M(M Q ) < 10 15 , 
with a relative mass distribution function given by (Q~|). The 
halos are assumed to be spatially distributed randomly with 
respect to each other. 

- C. Same as B, except that a fraction / su b of the matter in 
each halo of mass M is redistributed in the form of NFW 
subhalos in the mass range M x 10~ 4 < M sub < M x / sub , 
following relations ([8]), (0 and ( TTTT ). 

- D. Same as C, except that the parent halos are assumed to 
follow the halo mass function (Q]i in an absolute sense, so 
that only a fraction /haio(z) < 1 is in the form of halos at 
each redshift, whereas the remaining fraction at each red- 
shift is assumed to be smoothly distributed / sm ooth(z) = 1 - 
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Fig. 1. Schematic illustrations of matter distribution scenarios A-F. Here, the cylinder represents the cosmological volume be- 
tween the observer and source. In scenario A, all the matter is smoothly distributed along the sightline, in accordance with the 
Press & Gunn ( 1973 ) approximation. In scenario B, all the matter is instead assumed to be locked up inside halos with random 
spatial distribution. The effects of hierarchical halo clustering is here illustrated as evolution in the typical halo sizes between 
the source (high redshift; only low-mass halos) and observer (low redshift; both low- and high-mass halos). In C, substructure 
has been added to each halo. In D, a component of smoothly distributed matter (light gray) has furthermore been included to 
represent the matter that does not seem to be associated with halos within the mass range considered. Scenario E is similar to D, 
except that a halo mass catalogue from an N-body simulation has been used. Because of the resolution limit of this simulation, 
the number of small-mass halos has decreased. Scenario F takes into account the spatial clustering of halos predicted by N-body 
simulations, but is otherwise identical to scenario E. 



/haio(z). With a halo mass range of 10 10 < M(M Q ) < 10 15 , 
/smooth(z) increases from 0.55 at z ~ 0.2 to 0.96 at z ~ 5. 

- E. Same as D, except that the parent halos obey the mass 
function given by halo catalogues from the GIF N-body 
simulation (Kauffmann et al. 1999). Because of resolution 
effects, this limits the mass range to 10" < M(M Q ) < 10 14 . 
In this case, / smoo th(z) increases from 0.57 at z » 0.2 to 0.97 
atz ~ 5. 

- F. Same as D, except that the parent halos both obey the 
mass function and the spatial halo distribution given by halo 
catalogues from the GIF N-body simulation. 

The scenarios A-F are schematically illustrated in Fig. Q] 
Some of these levels of refinements correspond to assumptions 
adopted in the previous lensing literature. Scenario B is similar 
to the matter distribution model used in the lensing simulations 
by Bergstrom et al. (2000) for the magnification distribution 
of supernovae type la, whereas D is similar to the scenarios 
adopted by Oguri (2006) to study the image separation distri- 
bution in macrolensed sources and by Chen et al. (2003 to 
study the optical depth of small-mass halos along the line of 
sight to multiply-imaged quasars. 

The matter distribution adopted by Wyithe & Turner 
(2002a) in their investigation of the shortcomings of the Press 
& Gunn approximation (scenario A) when estimating the dis- 
tribution of microlensing optical depths, does not really corre- 
spond to any of the scenarios B-F, as they assume that no more 
than a single dark halo is located along each line of sight. We 
argue that this approach has severe limitations. In Fig. [2] we 
plot the average number of halos and subhalos (with scatter in- 



dicated by errorbars) located along the line of sight in scenario 
D, as derived from 1000 sightline simulations. 

Even for sources at redshifts as low as z s = 0.25, there are 
typically w 7 halos and w 2 subhalos along each sightline. As 
the source redshift is increased, these numbers grow monotoni- 
cally, reaching sa 40 for halos and « 10 for subhalos at z s = 5.0. 
It should be noted, however, that these estimates include every 
halo for which the sightline passes inside r v i r . As the sightline 
will traverse most of these halos at large distances from the cen- 
tre, where the density is low, the majority of halos will actually 
give a very small contribution to the surface mass density of the 
sightline. When sightlines do happen to pass through the dens- 
est inner regions of a halo, the contribution from that object to 
the surface mass density can on the other hand become con- 
siderable. Hence, it is in principle still possible that the surface 
mass density of entire sightlines may be dominated by individ- 
ual halos. To investigate how common this is, we isolate the 
single most dominant halo along each sightline and calculate 
its relative contribution to the overall surface mass density. The 
resulting probability distribution for the surface mass contribu- 
tion of the dominant halos is plotted in Fig.[3]for scenario D. 
The resulting probabilities for single-halo dominance are small 
even for nearby light sources, and drop rapidly with increasing 
source redshift. At z s = 0.25, the probability of having a sin- 
gle halo contributing more than 50% to the surface mass den- 
sity is about 6%, whereas already at z s = 1.0 this probability 
has dropped below 2%. In scenario D, these low probabilities 
arise partly because of the assumption of a non-zero fraction 
of smoothly distributed matter along each sightline, but even in 
scenario C, which assumes that all matter is locked up in halos 
(i.e. /smooth (z) = 0), the corresponding probabilities for single- 
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Fig. 2. The average number of halos (solid) and subhalos 
(dashed) along random sightlines to sources at redshifts z s in 
mass distribution scenario D. The error bars indicate the 68% 
probability limits on the sightline-to-sightline scatter around 
the average numbers of halos and subhalos, respectively. 

halo dominance are 26% at z s = 0.25 and 5% at z s = 1.0. 
Hence, individual halos are unlikely to dominate the surface 
mass density of sightlines to light sources at cosmological dis- 
tances. 

It must of course be emphasized that this conclusion holds 
only for random sightlines. In the study of macrolensed quasars 
- i.e. the very rare (roughly 1 in 500) objects that happen to 
have a foreground halo almost perfectly aligned along the line 
of sight - the situation can be quite different, as demonstrated 
by the investigation by Wambsganss et al. (2005), who find that 
90% of the sightlines to strongly lensed sources at z s ~ 1 are 
dominated by mass concentrations at a single redshift. 

3. Optical depth 

The microlensing optical depth t towards a light source at red- 
shift z s is defined as: 



At 



t = I cr(M, z, z s ) n(M, z)c — dz, 
Jo dz 



(15) 



where cr(M, z, z s ) is the cross section for a microlens of mass M 
at redshift z, and n(M, z) is the corresponding number density 
of such objects. Formally, the microlensing optical depth repre- 
sents the average number of lenses along a random line of sight. 
Under the assumption that the lenses do not overlap along the 
sightline, it also represents the fraction of sky that is covered 
by regions in which a point source will be microlensed. In the 
limit of small r - as e.g. in the local Universe - the optical depth 
can therefore directly be used as an estimate of the microlens- 
ing probability. At higher t, this interpretation does however 
break down because of overlapping microlenses. While the op- 
tical depth is difficult to directly relate to observations, it is of- 
ten used for quick-and-dirty estimates of the viability and rele- 
vance of different microlensing scenarios. Here, we will use it 



Q. 




single 



Fig. 3. Cumulative probability distribution p(< /single) for the 
relative surface mass density contribution /single from the sin- 
gle most dominant halo along a random sightline in matter dis- 
tribution scenario D. The different lines correspond to source 
redshifts of z s = 0.25 (thin solid with dots), 0.5 (thick solid), 
1.0 (thin solid), 2.0 (thick dashed), 3.0 (thin dashed), 4.0 (thick 
dash-dotted) and 5.0 (thin dash-dotted). 

to demonstrate the effect of different assumptions concerning 
the spatial distribution of MACHOs along the line of sight. 

In a cosmological model with flat geometry and £2 a + Om = 
1 .0, df/dz is given by: 

* = -1 l - . (16) 

dz #o(l+z) VQ m (1+z) 3 +O a 

If the microlenses are assumed to be isolated point-masses, the 
cross section for microlensing can be related to the Einstein 
radius Re through: 

AttGMD Xs D 0{ 



o-(M,z,Zs) = nR E 



c 2 D 



(17) 



where As, D \ and D os are the angular-size distances from lens 
to source, observer to lens and observer to source, respectively, 
for a lens at redshift z and a source at z s - Under the assump- 
tion that the MACHOs are uniformly and randomly distributed 
along the sightline (Press & Gunn 119731 scenario A), their 
number density n(M, z)pg is simply given by: 

£2macho(1 +z) 3 Pc 



«(M,z)pg 



M 



(18) 



By combining (fT5l>— (fT8T> we arrive at the optical depth r PG 
expected from this scenario: 

3// n M ACHO r Zs (1 +z) 2 D h D o] dz 



TpG 



(19) 



2D os Jo V q m(1 +z) 3 +O a 

By adopting filled-beam distances (see Kayser et al. 119971 
for a general discussion and Turner et al. 1984 for examples 
of how this assumption affects the lensing optical depth), the 
angular size distances are given by: 

dz 



(l+Zy)H J Zx ^Q m(1 +z )3 + Q a 



(20) 
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Fig. 4. The lensing optical depth expected from the Press & 
Gunn approximation (solid) compared to the distribution ex- 
pected from sightline statistics for different source redshifts, 
using matter distribution scenario B (all matter in halos; no 
subhalos). The errorbars indicates the standard deviation cr T of 
a log-normal distribution fitted to the predicted optical depth 
distribution. 



where z v and z v represent the redshifts of the points x and y. 

Because of the clustering of matter, the estimate < fT9b does 
not necessarily represent the true MACHO optical depth of ev- 
ery single sightline. In fact, most sightlines through the large- 
scale CDM structure of the Universe will appear underdense 
compared to the cosmic average. In situations where basically 
just a single line of sight is probed (e.g. when monitoring light 
sources located behind a particular galaxy or galaxy cluster), it 
may therefore be very useful to have a handle on the expected 
probability distribution of the optical depth. 

To numerically derive this distribution, we discretise the 
redshift axis between observer and source and compact the rel- 
evant intervening volume into i lens planes, so that the optical 
depth becomes: 



cal depth distribution at each source redshift can be reasonably 
well approximated by a log-normal distribution: 



t - ^ fzjo~(M, Zi, z s )n(M, z,)Z); 



(21) 



where fyj is given by ( fT4b . 

As expected, relaxing the Press & Gunn approximation 
causes substantial sightline to sightline scatter in t. We demon- 
strate this effect in Fig. |4]by comparing the optical depth Tpc 
derived under the Press & Gunn approximation to the distri- 
bution of optical depths resulting from a simulation of 1000 
sightlines at different source redshifts, using (f2Tb and mat- 
ter distribution scenario B. Unless stated otherwise, the cos- 
mological density of MACHOs is hereafter assumed to be 
^macho = = 0.3. Since the mass-dependence of n and <x 
cancel in (fl3T >. the assumed mass of an individual MACHOs is 
of no importance here. We find that, for this scenario, the opti- 



f(T,fi,cr iBT ) = 



1 



V27TTcr lnT 



exp 



(lnr-//) 2 



2cr. 2 
In? 



(22) 



with standard deviation cr\ aT and fi - lnrpc - cr 2 nr /2. As seen 
in Fig. |4] <r lnT decreases with increasing z s , The reason for this 
is the growing number of intervening halos at high redshifts, 
which reduces the relative scatter in r. While the Press & Gunn 
approximation clearly improves with source redshift, the scat- 
ter remains noticeable even at z s = 5. 

3.1. What is the most realistic matter distribution 
scenario? 

The predicted optical depth distribution at a given source red- 
shift is sensitive to the details of the matter distribution scenario 
assumed. While option F may seem like the most sophisticated 
approach out of the six scenarios described in Sect. 2.5, our 
current implementation of this scenario is limited by the rather 
high halo mass resolution limit (~ 10 11 M ) of the N-body 
simulation used. In Fig. [5] we compare the optical depth distri- 
bution predicted by scenario F to that of scenario E, in which 
identical assumptions concerning the mass function of dark ha- 
los are made, but in which the spatial distribution of the halos 
is random. As there is essentially no difference between the op- 
tical depth distributions predicted by these scenarios at any of 
the source redshifts investigated, we conclude that - at least 
down to the resolution limit of the N-body simulation used - 
the effects of halo clustering on the microlensing optical depth 
distribution is very small and can be neglected. The vertical part 
of the cumulative optical depth distributions seen at z s = 0.25 
stem from the fact that, with the limited mass resolution of 
the N-body simulations used, there is a non-negligible prob- 
ability of having zero halos along the sightline to low -redshift 
sources. For these sightlines, the optical depth is instead com- 
pletely dominated by the uniformly distributed matter fraction. 
This effect is also responsible for the artificial absence of tails 
towards small optical depths at higher redshifts and prevents 
the resulting distribution from being well-described by the log- 
normal fitting function d22"l >. 

In Fig. [6] we compare the microlensing optical depth dis- 
tribution predicted by scenarios A-D for a source at redshift 
z s = 1.0. As expected, the largest deviations from the Press 
& Gunn prediction (A) are produced by scenarios B and C, 
which assume all matter to be clustered into 10 I0 -10 15 M Q ha- 
los (i.e. /smoothCz) = 0). These two scenarios, which differ only 
in the inclusion of subhalos in C, actually produce very simi- 
lar optical depth predictions, indicating that subhalos represent 
a superfluous ingredient in high-redshift microlensing models. 
This result holds throughout the entire range of source redshifts 
investigated here (z s = 0.25-5.0). 

Since both halo clustering (in scenario F) and subhalos (in 
C) seem to have a negligible impact on the results, and since 
the N-body simulations used have limited halo mass resolution 
(E), the most useful computational scheme appears to be either 
scenario B or D. While B formally overestimates the cluster- 
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Fig. 5. The cumulative probabilities of having a sightline with 
microlensing optical depth lower than r in mass distribution 
scenarios E (random spatial distribution of halos; dashed) and 
F (spatial distribution of halos predicted by the N-body sim- 
ulation; thin solid) for source redshifts of (from left to right) 
z s — 0.25, 0.5 and 1.0. The vertical, thick solid lines indicate 
the corresponding optical depths predicted by scenario A (ran- 
domly and uniformly distributed compact objects). The negli- 
gible difference between the predictions of scenarios E and F 
indicate that, at least down to the resolution limit of the N-body 
simulation used, the effects of halo clustering on the microlens- 
ing optical depth distribution are very small. 



ing of matter into high-mass halos, scenario D, which includes 
a non-zero fraction of matter distributed uniformly, involves a 
certain degree of artificial smoothening. Because of these com- 
putational limitations, neither scenario appears perfectly sat- 
isfying in generating realistic predictions for the line-of-sight 
density distributions expected in a ACDM Universe. In what 
follows, we therefore explore the consequences of both scenar- 
ios B and D, as the most realistic situation should lie some- 
where in between. 

The optical depth distribution derived for our matter dis- 
tribution scenarios are quite different from those presented by 
Wyithe & Turner (2002a). Since most sightlines are under- 
dense, the median r of our distributions is always lower than 
tpg. The skewness of the distributions derived by us is, how- 
ever, not nearly as extreme as that predicted by Wyithe & 
Turner. In their model, the vast majority (90-98%) of all sight- 
lines have t < Tpc, whereas in our case, the corresponding frac- 
tion is only 60-80%. This discrepancy can be traced to their as- 
sumption of having only a single halo located along each sight- 
line. As demonstrated in Figs.[2]and[3] dark halos overlap sub- 
stantially, and since their cross sections are dominated by their 
extended, low-density outskirts (i.e. regions with small t), the 
dominating effect of including more than one halo along each 
sightline is to reduce the probability for very small r, thereby 
shifting the t distribution in the direction of the tpg estimate. 




Fig. 6. The cumulative probabilities of having a sightline with a 
microlensing optical depth lower than r for a source at z% = 1.0 
in matter distribution scenarios B (all matter in halos and no 
subhalos; thin solid), C (all matter in halos and subhalos; thick 
dashed) and D (a fraction of matter in halos and subhalos; thin 
dashed). The vertical, thick solid lines indicate the optical depth 
predicted by scenario A (randomly and uniformly distributed 
compact objects). 

3.2. The predicted distribution of MACHO optical 
depths 

The redshift dependence of the log-normal standard deviations 
cn nT derived from the optical depth distribution predicted by 
scenarios B and D are well-described by the following polyno- 
mials: 

o-i nT>B (z s ) = -0.004015^ + 0.05854z* - 0.3232z^ (23) 
+0.8539z^- 1.156z s + 1.014; 



o-inr.Dfc) = 0.00041994 - 0.00090954 - 0.02956z^ (24) 
+0.1833z^ - 0.4042z s + 0.4993. 

Using the fitting functions (23) and (24), the approxi- 
mate optical depth distribution around the analytical estimate 
( |T9t can easily be calculated, thereby allowing the viability 
of the Press & Gunn approximation to be assessed in various 
microlensing situations involving cosmologically distributed 
MACHOs. Please note that the fitting functions are only valid 
for source redshifts in the range z s = 0.25-5. Because of the 
negligible impact of subhalos, scenario C produces a fitting 
function that is virtually identical to that of D. Due of the lim- 
ited halo mass resolution in scenarios E and F, the optical depth 
distributions generated by these are not very realistic, and can- 
not be well-reproduced by a lognormal function. 

In scenario D, one standard deviation corresponds to a vari- 
ation on the optical depth of « 50% for low-redshift (z s = 0.25) 
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light sources. While the scatter drops substantially with in- 
creasing source redshifts, it remains noticeable (at the w 10% 
level) even at z s = 5.0. In scenario B, the corresponding scatter 
is even higher. At z s = 0.25, the standard deviation corresponds 
to a factor of « 2.2, and even at z s = 5.0, the l<x scatter is at the 
« 25% level. 

Because of the redshift dependence of the dark halo mass 
function, the fraction of matter that is within the halo mass res- 
olution limits of scenario D decreases as a function of redshift. 
Whereas about half of the matter content is locked up in re- 
solved halos at z = 0, the corresponding fraction at z — 5 is 
only a few percent (see Sect. 2.5). The difference between the 
matter distributions of scenarios B and D therefore becomes 
more pronounced with increasing z s . This does not, however, 
result in any dramatic difference in the z s dependence of cr\ nT 
predicted by the two scenarios, since the relative contribution 
to the optical depth from matter at the highest redshifts is small. 
For example, a light source at z s =5.0 receives only 16% of its 
optical depth from matter at z > 3 and only 4% from matter at 
z > 4. 

The error on the tpg estimate may become important when 
attempting to derive Qmacho from microlensing observations 
of only a few unique sightlines, as e.g. in the case of surface 
brightness fluctuations measured in galaxies selected from a 
small patch of the sky (Lewis & Ibata 2001). In the limiting 
case of only a single sightline, cr\ nT represents an uncertainty 
below which Qmacho cannot be determined without attempting 
to reconstruct the detailed surface mass density contribution of 
that particular sightline. A simple way to avoid this uncertainty 
is of course to select targets along widely separated, uncorre- 
lated sightlines. The observational programmes proposed by 
Tadros et al. (2001 ) and Totani (2003) to monitor sources be- 
hind galaxy clusters, do for instance not suffer from this prob- 
lem, since the low redshifts of their targeted foreground objects 
imply that sources will be selected over very large angles. 

The error on the tpg estimate may also affect, albeit to 
much smaller extent, estimates of the relative contribution from 
cosmologically distributed MACHOs to the microlensing vari- 
ability displayed by a single sightline crossing a dominant mass 
condensation with moderately high microlensing optical depth. 
As an example, consider the case of a light source at z s =0.5 lo- 
cated behind a foreground halo (galaxy or galaxy cluster) with 
an impact parameter such that the MACHO optical depth of 
that halo is expected to be r ha i = 0.2 /macho, where /macho 
represents the fractional contribution of MACHOs to the matter 
content of the halo. At z s = 0.5, the optical depth of cosmolog- 
ically distributed MACHOs is, in the Press & Gunn scenario, 
tmacho ~ 0.06 Qmacho (neglecting magnification bias) or 
~ 0.018 /macho (for Q M = 0.3), if we assume that the MACHO 
mass fraction /macho is scale-independent. Hence, the op- 
tical depth contribution from the cosmologically distributed 
MACHO population is expected to be one order of magnitude 
smaller than that of the MACHOs in the foreground object. If 
the total microlensing optical depth is observationally deter- 
mined to be r bs = 0.08 (with negligible observational error), 
one would conclude that /macho = 0.08/(0.2 + 0.018) ~ 0.37. 
The uncertainty cr\ nT » 0.6 predicted for this source redshift 
in scenario B does however translate into a 95% confidence in- 



terval of /macho = 0.37jwj? on this measurement. While this 
error is admittedly small, it will be compounded by the error 
on the Thaio estimate and may become relevant in future high- 
precision observations aimed to constrain /macho- 

4. Application to the long-term optical variability 
of quasars 

While the optical depth is a concept that is easy to understand 
and useful for demonstrating the effect of various mass distri- 
bution scenarios, it is somewhat difficult to directly relate to 
observations. One reason for this is that observable microlens- 
ing signatures tend to depend not only on the number of lenses 
affecting a sightline, but also on their velocities and where 
along the sightline the lenses are located. Another reason is 
that the optical depth, as defined here, represents the number of 
lenses along a random sightline, whereas flux-limited observa- 
tions probe magnification-biased (rather than random) lines of 
sight. The calculations presented so far furthermore only take 
into account the effects of microlenses with a projected dis- 
tance below one Einstein radius from the line of sight, whereas 
non-negligible magnification may in fact also result from the 
accumulated effects of lenses with larger impact parameters. 

To demonstrate how the interpretation of observational data 
is affected by the use of the Press & Gunn approximation, 
we couple the matter distribution model described in Sect. 2 
to a microlensing code and revisit the long-standing issue of 
whether the optical variability of quasars on long time scales 
(years to decades) can be caused by microlensing by MACHOs 
in the subsolar mass range. In principle, this microlensing vari- 
ability scenario has several attractive features. It provides a nat- 
ural explanation for the statistical symmetry (Hawkins 2002), 
achromaticity (Hawkins 2003) and lack of cosmological time 
dilation (Hawkins 2001) in the optical light curves of quasars. 
In Zackrisson et al. (120031 1 we did however demonstrate that, 
contrary to the claims by Hawkins, the optical, long-term vari- 
ability of quasars at redshifts z s < 1 cannot be attributed to 
microlensing alone. Some other variability mechanism - most 
likely accretion disk instabilities - must also contribute sub- 
stantially. While this does not necessarily rule out the possibil- 
ity that microlensing may contribute at some level, it compli- 
cates the procedure of using light curve statistics to infer the 
properties of any potential MACHO populations along the line 
of sight. 

4.1. Microlensing model 

To simulate the lensing situation relevant for the long-term op- 
tical variability of quasars, we couple our matter distribution 
model to the microlensing model described in Zackrisson & 
Bergvall (2003). In short, this microlensing code distributes 
compact objects in a number of lens planes between the source 
and observer. After each time step, the lenses are repositioned 
according to their individual velocities. The individual magnifi- 
cation contributions from each separate microlens is computed 
using analytical approximations and their total magnification 
computed using the multiplicative magnification approxima- 
tion (Ostriker & Vietri 1983). The most important model pa- 
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Fig. 7. The cumulative probabilities P(< 6m) of observing 
quasars with amplitudes lower than 6m predicted in the case 
of z s = 0.5, Qmacho = 0.3, Rs = 10 12 m and lens masses (from 
top to bottom at 6m = 0.3) Mmacho = 10~ 2 M Q , 10~ 3 M Q and 
10 M Q . For each lens mass, the different lines correspond 
to the results from matter distribution scenarios A (solid), B 
(dashed) and D (dash-dotted). 

rameters are: the parameters describing the background cos- 
mology (O m , Q A , Hq and the homogeneity parameter rj); the 
parameters describing the compact object population (Qmacho 
M M acho an d their velocity dispersion <x v macho perpendicu- 
lar to the line of sight) and parameters describing the source 
(source radius and redshift zs), which is assumed to be a 
disk of uniform brightness. The limitations of this microlens- 
ing code are thoroughly described in Sect. 5 of Zackrisson & 
Bergvall d2003l) and Sect. 7 of Zackrisson et al. (TZOOH) . 

The matter distribution scenarios described in Sect. 2. are 
implemented in the Zackrisson & Bergvall (2003) by multiply- 
ing the density of microlenses p, ; macho,pg in each lens plane i, 
expected in the framework of the Press & Gunn approximation 
(scenario A), by the factor fej (changing stochastically from 
sightline to sightline) defined in (TBl i: 

P/.MACHO = /s,/P/,macho,pg (25) 

While the lensplane-to-lensplane scatter in number of mi- 
crolenses changes dramatically with this approach, we still as- 
sume a random spatial distribution of the MACHOs located 
in each lens plane. Hence, we assume that there is no ad- 
ditional small-scale clustering on scales below the resolution 
limit of our matter distribution model (e.g. in the form of bi- 
nary MACHOs or MACHOs bound into small clusters). 

4.2. The amplitude distribution and its redshift 
dependence 

To assess the contribution of MACHO microlensing to the op- 
tical long-term variability of quasars, the relevant quantity to 
predict is the level of brightness variability that a particular 



compact object population will induce during the course of an 
observational monitoring programme. Here, we will concen- 
trate on the light curve amplitudes, Sm, defined as the difference 
between the minimum and maximum of the yearly mean mag- 
nitudes m(t) observed when a source is monitored for a certain 
period of time with a particular sampling rate: 

6m = max(m(f)) - min(m(f))- (26) 

In Hawkins (f2000), this quantity was derived for several well- 
defined samples of quasars that had been monitored several 
times a year for 20 years. As demonstrated in Zackrisson et 
al. (2003), microlensing predicts too small mean amplitudes 
at z s < 1 and an overall amplitude distribution with too few 
6m > 0.35 mag quasars. 

To predict the amplitude statistics of microlensing-induced 
variability, we generate a large number of light curves with a 
time sampling that closely matches those of the observational 
samples of Hawkins d20001 l. To each yearly magnification data 
point, derived by taking the average of four evenly spaced in- 
trayear magnification data points, we then add Gaussian noise 
with standard deviation <x m = 0.065, which represents a com- 
bination of observational error and short-term variability (as- 
sumed to be intrinsic). From these yearly magnification data 
points, the amplitude is then derived. 

As the observational samples are flux limited, we further- 
more assume a quasar luminosity function and assign each light 
curve to a B-band absolute magnitude Mb - This quantity is con- 
verted into an apparent magnitude m,B using the magnification 
of the first data point of the light curve and an assumed k- 
correction. All light curves that fail to meet the magnitude limit 
of the observational sample are then rejected. Because there is 
a correlation between the light curve amplitude and the mag- 
nification used to derive the apparent magnitude of each syn- 
thetic quasar (Schneider |19931 l, this procedure imposes a slight 
magnification bias. In the following, the Mb = 21.5 magni- 
tude limit of the Hawkins (2000) UVX sample will be assumed. 
Zackrisson & Bergvall (2003) and Zackrisson et al. (2003) can 
be consulted for further details on the sample simulation pro- 
cedure. 

As noted by Zackrisson et al. d20031 l. one naively expects 
that relaxing the Press & Gunn approximation would widen 
the distribution of light curve amplitudes, as quasars located 
along overdense sightlines should display increased variabil- 
ity, whereas quasars located along underdense sightlines should 
vary less. It is not obvious, however, that the effect will be very 
pronounced, since the amplitudes are dominated by the masses, 
velocities and impact parameters of the microlenses - quanti- 
ties which are unaffected by the matter distribution model ex- 
plored here. What is affected is instead the probability distribu- 
tion for having several microlenses simultaneously contribut- 
ing substantially to a light curve. 

We find, that while relaxing the Press & Gunn approxima- 
tion does indeed affect the amplitude distribution in the way 
foreseen by Zackrisson et al. (2003), the effect is typically 
small and varies in a non-trivial way across the relevant mi- 
crolensing parameter space. The main reason for this is that 
when the light curve is characterized by distinct microlens- 
ing peaks, rather than multiple microlensing peaks superim- 
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Fig. 8. The cumulative probabilities P(< 6m) of observing 
quasars with amplitudes lower than 6m predicted for Qmacho = 
0.3, R s = 10 12 m, M = 1(T 4 M@, cr v = 400 km s _1 in the case 
of matter distribution scenarios A (solid), B (dashed) and D 
(dash-dotted). From left to right, the three different sets of lines 
correspond to source redshifts of z s = 0.25, 0.5 and 1.0. 



posed on each other, the amplitude distribution is completely 
insensitive to variations in the line-of-sight density around the 
mean. As long as the density variations are not large enough to 
substantially alter the probability for peak superpositions, the 
Press & Gunn approach (scenario A) therefore gives a perfectly 
adequate description of the amplitude distribution. When, on 
the other hand, the amplitudes of essentially all sightlines are 
dominated by overlapping microlensing peaks, density varia- 
tions do have a slight effect on the amplitude probability dis- 
tribution. The transition between these two regimes is illus- 
trated in Fig. [7] where we plot the amplitude distribution for 
a light source at z s = 0.5 and MACHO masses in the range 
10 -4 — 10 2 M Q , under the assumption of Qmacho = 0.3 and 
Rs = 10 12 m. For a M M acho > 10~ 2 M Q object, the light 
curves typically only contain a single dominant microlensing 
peak, and the amplitude predictions of scenarios B and D are 
indistinguishable from those of scenario A. Hence, the Press 
& Gunn approximation can safely be used to model the effects 
that the ~ 10 _1 M objects that may have been detected in the 
halos of the Milky Way and M31 (Alcock et al. 120001 Calchi 
Novati et al. 2005) would have on these amplitude distribu- 
tions. At 10~ 3 M Q , a small, but noticeable difference does oc- 
cur, and at 10~ 4 M , the lines are clearly separated, albeit not 
by much. While the average amplitudes are essentially iden- 
tical in the different scenarios, the amplitude distribution does 
become somewhat wider as a result of adopting a more realistic 
spatial distribution for the MACHOs. 

The source size assumed here {Rs = 10 12 m) is smaller than 
that typically used in the literature by an order of magnitude, 
and reflects the value found by Zackrisson et al. (12003) to be 
the most promising for explaining the high level of variability 
observed at low redshifts (z s < 1) in the UVX quasar sample of 



Hawkins (2000). This choice allows us to assess the prospects 
of microlensing as an explanation for the observed variability, 
in light of the more realistic matter distibution scenarios dis- 
cussed in the current paper. Our tests indicate that going from 
Rs = 10 12 m to the canonical value Rs = 10 13 m only has a 
modest impact on the relative difference between the amplitude 
distributions generated by scenarios A, B and D at the MACHO 
masses explored in Fig. [7] 

If the microlensing parameters cr v macho Qmacho and Rs 
are kept fixed, the balance between single- and multiple-peak 
microlensing moreover depends on source redshift, as lower 
z s correspond to a lower optical depth. This is illustrated in 
Fig [8] where the amplitude distribution predicted in the case of 
a Mmacho = 10~ 4 M Q , Qmacho = 0.3, R s = 10 12 m scenario 
is plotted for source redshifts z s = 0.25, 0.5 and 1.0. This cor- 
responds to the microlensing scenario found to be one of the 
most promising for explaining the long-term optical variability 
of quasars (Zackrisson et al. 2003 1 but note that no microlens- 
ing scenario was found to be completely successful in this re- 
spect). This scenario is admittedly something of a stretch, as 
it implies that all of the matter in the Universe (baryonic and 
non-baryonic) is in the form of MACHOs. This particular pa- 
rameter combination has moreover been ruled out by the up- 
per limits on the cosmological density of compact objects de- 
rived by Zackrisson & Bergvall (2003) and Zackrisson et al. 
(2004). Since these investigations employed the Press & Gunn 
approximation to reach these conclusions, it is nonetheless im- 
portant to investigate how the conclusions from investigations 
are affected when this approximation is relaxed. The MACHO 
mass assumed here is intermediate between that advocated by 
Schild ([19961 Pelt et al. (fT998l and Colley & Schild (120031 
for explaining the rapid variability in Q0957+561, and the es- 
timate advocated by Hawkins (1996) for explaining the long- 
term variability of quasars that are not multiply-imaged. While 
certainly not favoured by the constraints on MACHOs imposed 
by the MACHO and EROS/EROS2 surveys (e.g. Lasserre et al. 
2000; Tisserand et al. 2007), diffuse objects of this mass may 
still be viable, as discussed in the introduction. At z s = 0.25, 
the difference between scenarios A and B/D is very small, even 
though the density variations along the line of sight are sub- 
stantial, simply because the number of lenses affecting the light 
curve is not large enough for superimposed microlensing peaks 
to dominate the amplitude distribution. As shown in Fig. [8] 
the difference between the predictions of the different scenar- 
ios grows slightly with increasing z s , as an increasing fraction 
of the sightlines become dominated by multiple events, even 
though the relative sightline density scatter goes down. 

In Fig. [9] we display the amplitude distribution predicted 
by scenarios A (thin solid), B (thin dashed) and D (thin dash- 
dotted) for a distribution of quasar redshifts that matches that 
of the UVX sample from Hawkins (2000), and contrast the pre- 
dictions to the observed amplitude distribution of this sample 
(thick solid). While the effect is small, the trend is - as sug- 
gested by Zackrisson et al. (2003) - to increase the discrep- 
ancy between microlensing simulations and the observations, 
which is most notably seen by the lack of microlensing light 
curves with amplitudes 6m > 0.35 magnitudes. Hence, the 
problems associated with microlensing as the dominant mech- 
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Fig. 9. The cumulative probabilities P(< 6m) for amplitudes 
lower than 5m in the UVX sample of Hawkins d2000t thick 
solid), compared to the predictions for Qmacho = 0.3, Rs = 
10 12 m, M = 1(T 4 M , cr v = 400 km s -1 and scenario A (thin 
solid), B (thin dashed) and D (thin dash-dotted). Here, a source 
redshift distribution matching that of the UVX sample has been 
used. While the effect of relaxing the Press & Gunn approxi- 
mation (i.e. going from scenario A to B or D) is modest, it 
does increase the fraction of amplitudes 5m < 0.35 mag by a 
few percent, and therefore increases the discrepancy between 
observations and microlensing predictions. 

anism for the long-term variability are in no way diminished, 
but instead slightly augmented, once the large-scale clustering 
of MACHOs is taken into account. The implementation of ad- 
ditional features of MACHO clustering (as described in Sect. 
5) are expected to drive the amplitude distribution predicted by 
microlensing even further away from the observed one. Since 
the effect of relaxing the Press & Gunn approximation has been 
shown to be quite modest (on the level of a few percent), the 
Press & Gunn approach may however still be used to gener- 
ate MACHO microlensing predictions of reasonable accuracy 
for the amplitude distribution of quasars. This greatly reduces 
the computational complexity of MACHO microlensing mod- 
els for high-redshift sources. 

Throughout this section, we have for clarity plotted the pre- 
dictions of scenarios A, B and D only. As expected from the 
analysis of optical depth distributions in Sect. 3, the light curve 
amplitude predictions of scenario C are essentially identical to 
those of scenario B, whereas the amplitude distributions gener- 
ated by scenario E are indistinguishable from those of scenario 
F and fall between those of scenarios A and D. 

5. Discussion 

As we have shown, the Press & Gunn (U973I I approximation 
neglects the source-to-source scatter in microlensing properties 
for light sources at cosmological distances. While the approx- 
imation improves with increasing source redshift, the scatter 
remains detectable up to a redshift of at least z = 5, in the 



case where the microlenses follow the large-scale distribution 
of CDM. If part of the dark matter of the Universe truly is made 
up by MACHOs in the stellar mass range, as suggested by local 
microlensing detections (e.g. Alcock et al. 2000; Calchi Novati 
et al. 2005), microlensing by such objects is expected to pro- 
duce observable effects in a non-negligible fraction of high- 
redshift light sources like supernovae type la and gamma-ray 
bursts (e.g. Minty et al. 120021 Wyithe & Turner l2002bl Baltz 
& Hui 2005). Models of the type presented here may there- 
fore prove to be important tools in the statistical analysis of 
microlensing events detected in ongoing and upcoming obser- 
vations (with e.g. Swift, SNAP, JEDI or Destiny) of such ob- 
jects. 

The matter distribution model used here is more sophisti- 
cated than the one used by Wyithe & Turner (2002a) in that it 
takes into account the combined effect of several halos along 
the line of sight, their spatial clustering and internal substruc- 
ture. There is, however, still room for improvement, as there are 
a number of aspects of CDM clustering that our present model 
does not treat. 

The model assumes, to begin with, that all matter not 
locked up inside halos of mass larger than ~ 10 10 M (scenar- 
ios B-D) or ~ 10 11 M (scenarios E-F) is uniformly distributed 
along the line of sight. This neglects a) the spatial clustering 
of low-mass field halos which current cosmological N-body 
simulations have problems in tracing due to limited particle- 
resolution, and b) loosely bound matter in filaments or else- 
where that are left over by the halo finders used to create halo 
catalogues from N-body simulations. As this can correspond to 
a substantial fraction of the matter (> 50% in the case of the 
GIF simulation used here), this component should be modelled 
with more care. While the lower mass limit of the halo cata- 
logues can be lowered with increasing particle resolution, tak- 
ing filamentary, non-halo material into account will most likely 
imply using the raw N-body simulation data instead of halo 
catalogues and analytical density profiles, which will make the 
model substantially more computationally demanding. 

The current matter distribution model assumes that subha- 
los make up a fixed fraction of the material of each halo of 
a given mass, whereas the actual subhalo mass fraction dis- 
plays substantial scatter (see e.g. van den Bosch et al. 2005; 
Shaw et al. 120071 Nurmi et al. 120061) . As the very modest dif- 
ference between the results generated by our matter distribu- 
tions B and C already indicates that subhalos are unimportant 
for high-redshift microlensing statistics, this is not expected to 
have any dramatic effects on our results. 

By using the spherically averaged NFW density profiles, 
the effects of halo triaxiality and the associated scatter in sur- 
face mass density coming from different halo viewing angles 
are neglected. This situation can be improved by studying the 
distribution of surface mass densities produced when the line of 
sight traverses realistic N-body halos at a given impact parame- 
ter. If an analytical function can be fitted to this distribution, the 
increased scatter introduced by both triaxiality and substructure 
can then be statistically reproduced by Monte Carlo techniques 
for each NFW sphere along the sightline in the matter distribu- 
tion model. We are currently investigating the feasibility of this 
approach (Holopainen et al. 2007). 
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Our matter distribution model is indirectly based on results 
from dissipationless N-body simulations. In reality, however, 
CDM halos contain baryons, which are expected to alter both 
the shapes and density profiles of the dark halos. The general 
expectation is that baryon cooling inside CDM halos will cause 
the dark halos to contract (Blumenthal et al. [1986; Gnedin et 
al. 120041 Sellwood & McGaugh [2005b . although bar forma- 
tion could possibly drive the evolution of the innermost part of 
dark halo in the opposite direction (e.g. Holley-Bockelmann 
et al. 120051 but see Colin et al. 120061 for a different view). 
Baryon cooling is furthermore expected to render dark halos 
more spherical (e.g. Kazantzidis et al. 2004 ), but the exact mag- 
nitude of this effect remains controversial. 

At the current time, the velocities of the microlenses are as- 
sumed to follow a Gaussian distribution which is independent 
of their spatial location. In principle, this approximation can be 
improved by using the analytical or numerical density profile 
of each halo to calculate the velocity distribution of compact 
objects in the relevant volume of that object, and then convolve 
this distribution with the peculiar velocity of each halo, which 
may be derived from N-body simulations. For matter not as- 
sociated with halos, the problems becomes substantially more 
complicated. 

The current model also neglects spatial clustering of mi- 
crolenses on scales below the smallest structures considered 
(the least massive subhalos of mass 10 6 M Q ). If MACHOs 
are bound into binary systems or small clusters, this will af- 
fect both their position and velocity distributions in each lens 
plane. Without a detailed prediction for the small-scale distri- 
bution of MACHOs, it is however difficult to fully quantify the 
magnitude of this effect. Very dense, small-scale systems of 
MACHOs are of course subject to a number of dynamical con- 
straints, as reviewed by Carr & Sakellariadou ( 1999 ). 

The microlensing model used here to assess the impact of 
various assumptions concerning the spatial distribution of com- 
pact objects along the line of sight is based on the multiplicative 
magnification approximation (Ostriker & Vietri [19831) . While 
this is not expected to challenge the overall results of the paper, 
the possibility that spurious effects may arise in certain parts 
of the microlensing parameter space (especially in the small 
source regime) as a result of this, cannot be ruled out, as dis- 
cussed at some length in Zackrisson & Bergvall (2003 ) and 
Zackrisson et al. {2003 ). This concern can be resolved by em- 
ploying ray-shooting techniques, at the expense of a substan- 
tial increase in required CPU-time. The current lensing cal- 
culations furthermore neglect contributions to the magnifica- 
tion bias from structures larger than the microlenses themselves 
(Pei ll995l ). As discussed in Zackrisson & Bergvall (120031) . the 
microlensing model used furthermore assumes filled-beam an- 
gular size distances, which is admittedly not a self-consistent 
treatment of the situation, and may be a bad approximation 
in certain parts of the microlensing parameter space. The 
sightline-to-sightline variations predicted by the matter distri- 
bution scenarios considered here act to augment this problem 
even further. To remedy this problem would however require 
a much more sophisticated lensing algorithm, e.g. of the Holz 
& Wald (1998} type, in which "each photon creates its own 
Universe" along its path from source to observer. For a fixed 



Qmacho. the effect of decreasing the filling factor of each beam 
(i.e. going from filled to empty) is to decrease the the opti- 
cal depth (e.g. Turner et al 1984). Because of this, the sight- 
lines that correspond to low microlensing optical depths in our 
formalism would tend to have them lowered even further if 
the filled-beam approximation was dropped, thereby increas- 
ing the sightline to sightline scatter. This effect is however not 
expected to be dramatic unless the source redshift is very high, 
as it is only then that the filled-beam and empty-beam angu- 
lar size distances differ substantially. By evaluating equation 
d"L9l when shifting between filled-beam and empty-beam angu- 
lar size distances, we find that the difference in optical depth 
between these two extreme situations (the most realistic situ- 
ation should lie somewhere in between) is < 3 % for source 
redshifts of z s < 1, which is negligible compared to the scat- 
ter induced by other components of the the matter distribution 
model. However, at z s =5 the difference is < 23 %, which 
is sufficiently large to warrant consideration in future, high- 
precision predictions of the optical depth distribution. 

6. Summary 

When predicting the gravitational microlensing effects of dark 
matter in the form of compact objects (MACHOs) on high- 
redshift light sources, it is customary to adopt the simplify- 
ing assumption that MACHOs are randomly and uniformly dis- 
tributed along the line of sight (the Press & Gunn ll973l approx- 
imation). In this paper, we explore the consequences of relax- 
ing this assumption in favour of a model where the microlenses 
follow the spatial clustering of CDM, as would be expected for 
scenarios involving non-baryonic MACHOs. In summary, we 
find that: 

- The contribution from cosmologically distributed 
MACHOs to the microlensing optical depth of indi- 
vidual sightlines may differ from the Press & Gunn 
estimate (T% by a non-negligible factor, which depends on 
the redshift of the light source studied. Here, we quantify 
the scatter as a function of source redshift and demonstrate 
that the optical depth of an individual sightline can be in- 
correct by a factor of * 2 (at the lcr level) for a light source 
at z s = 0.25. While the optical depth scatter decreases with 
source redshift, it remains noticeable even at z s = 5. This 
result implies that attempts to determine the cosmological 
density of MACHOs from microlensing observations of a 
single or only a few random sightlines will be subject to 
substantial uncertainties. 

- The sightline-to-sightline scatter of microlensing optical 
depth is dominated by the clustering of MACHOs into 
galaxy-mass CDM halos. Both the spatial correlation of 
these halos and their internal dumpiness due to substruc- 
ture is shown to have negligible impact on the optical depth 
statistics. 

- While relaxing the Press & Gunn ( 119731 ) approximation in 
favour of the matter distribution scenarios explored here is 
shown to have a rather modest impact on the MACHO mi- 
crolensing predictions for light curve amplitudes, the mi- 
crolensing predictions are driven further away from the ob- 
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served statistical properties of the long-term optical vari- 
ability of quasars that are not multiply-imaged. Hence, the 
problems with microlensing as the primary cause for this 
variability are in no way diminished, but instead slightly 
augmented, once the large-scale clustering of MACHOs is 
taken into account. 
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